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Abstract 

FeynHiggsFast is a Fortran code for the calculation of the masses and the mixing 
angle of the neutral CP-even Higgs bosons in the MSSM up to two-loop order. It is 
based on a compact analytical approximation formula of the complete diagrammatic 
one-loop and the dominant two-loop contributions. At the one-loop level a leading log- 
arithmic result is used, taking into account all sectors of the MSSM. At the two-loop 
level at 0{aa s ) the leading logarithmic and non-logarithmic contributions are taken 
into account. The approximation formula is valid for arbitrary choices of the parame- 
ters in the Higgs sector of the model. Comparing its quality to the full diagrammatic 
result, we find agreement better than 2 GeV for most parts of the MSSM parameter 
space. 
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1 Introduction 



One of the most striking phenomenological implications of Supersymmetry (SUSY) is the 
prediction of a relatively light Higgs boson, which is common to all supersymmetric models 
whose couplings remain in the perturbative regime up to a very high energy scale [0]. The 
search for the lightest Higgs boson thus allows a crucial test of SUSY and especially the 
Minimal Supersymmetric Standard Model (MSSM) [[J. Therefore it is one of the main goals 
at the present and the next generation of colliders. A precise knowledge of the dependence 
of the masses and mixing angles of the Higgs sector of the MSSM on the relevant SUSY 
parameters is necessary for a detailed analysis of SUSY phenomenology at LEP2, the up- 
graded Tevatron, and also for the LHC and a future linear e + e~ collider, where high-precision 
measurements in the Higgs sector might become possible. 

The mass of the lightest CP-even MSSM Higgs boson, m^, is bounded from above by 
about m h < 135 GeV, including radiative corrections up to the two-loop level [|]-|7j]. The 
mixing angle a e g that diagonalizes the neutral CP-even Higgs sector receives the same kind 
of corrections. By incorporating a e fr into the Higgs decay widths, the leading electroweak 
corrections to the decay of the neutral CP-even Higgs bosons are taken into account. 

In this note we present the Fortran code FeynHiggsFast. Using low energy MSSM pa- 
rameters as input, it evaluates the masses of the neutral CP-even Higgs bosons, nth and 
irtHi as well as the corresponding mixing angle, a c fr, at the two-loop level. In addition the 
mass of the charged Higgs boson, m#± , is evaluated at the one-loop level. The p -parameter, 
leading to constraints in the scalar fermion sector of the MSSM, is evaluated up to 0(aa s ), 
taking into account the gluon exchange contribution at the two-loop level ||. FeynHiggs- 
Fast is based on a compact analytical approximation formula, containing at the two-loop 
level the dominant corrections in 0(aa s ) obtained in the Feynman-diagrammatic (FD) ap- 
proach HU and subdominant corrections of 0[G 2 F m\ ) obtained with renormalization group 
(RG) methods PH7|. In comparison with the FD result, consisting of the complete one-loop 
and the two-loop contributions as given in Ref. [||] (incorporated into the Fortran code Feyn- 
Higgs ||), the approximation formula is much shorter. Thus FeynHiggsFast is about 3 x 10 4 
times faster than FeynHiggs. Agreement between these two codes of better than 2 GeV is 
found for most parts of the MSSM parameter space. 

The following sections are organized as follows: In Sect. 2 we shortly summarize the 
analytical approximation formula and perform a comparison with the full FD result. A 
description of how to use FeynHiggsFast is given in Sect. 3. The conclusions can be found 
in Sect. 4. 

2 The analytical approximation 

In order to fix our notations, we first list the conventions for the MSSM scalar top sector: 
the mass matrix in the basis of the current eigenstates ti and tn is given by (neglecting 
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Furtheron the following simplification is used: 



In this simplified case we defineQ 



M h =M- tR := m- q . 



M 2 s :=m 2 + m 2 t . 



For further details see Ref. 10 
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At the tree level the mass matrix of the neutral CP-even Higgs bosons in the 0i, 02 basis 
can be expressed as follows: 
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Ma denotes the mass of the CP-odd Higgs boson, Mz is the mass of the Z boson, and 
tan/? = V2/V1 is the ratio of the two vacuum expectation values of the two Higgs doublets 
in the MSSM, see Ref. 
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At higher orders, the Higgs mass matrix (|5]) is supplemented by the renormalized self- 
energies T, s (q 2 ) , s = 0i,02,0i02- Here we use the approximation of neglected external 
momentum: S s = E s (0). The masses and mn are then obtained by diagonalizing the 
higher order corrected mass matrix with the mixing angle a e ff: 
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Here we only give a very brief outline of the calculation of the renormalized Higgs boson 
self-energies. A detailed description can be found in Ref. [jlCfl . 

The mass of the lightest Higgs boson receives contributions from all sectors of the MSSM, 
but not all are numerically of equal relevance. The dominant corrections arise from the t — t 
sector of the MSSM: The results for the renormalized self-energies S s have been derived 
analytically in the FD approach in Refs. |3|,U. These results, however, are rather lengthy. 
In order to derive a compact analytical expression, several approximations have been made 



1 The more general case with Mf L ^ Mf R can be found in eq. (|T^ 
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as described in Ref. [10]. The main step of the approximation consists of a Taylor expansion 
in 

\m t X t \ , . 



A* 



M 2 



of the S s (0). For the one-loop correction the expansion has been performed up to 0(A~); all 
three renormalized one-loop Higgs-boson self-energies give a contribution. Concerning the 
two-loop self-energies, the expansion has been carried out up to 0(A~); in the approximation 

~ (2) 

considered here only S^ 2 gives a non-zero contribution. 

Leading contributions beyond 0(aa s ) have been taken into account by incorporating 
the leading two-loop Yukawa correction of 0(G 2 F m®) [§-0. Furthermore the result has been 
expressed in terms of the MS top-quark mass 



m t = m t lm t 
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1 + 7^a s (m t ) 



(9) 



instead of the pole mass m t . This leads to an additional contribution in 0{aa 2 s ). 

The analytical expressions arising from the t — t sector are given at the one-loop level as 
follows: 
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The two-loop contributions in 0(aa s ) read: 
Eg(0) = 0, 
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The two-loop Yukawa correction in this approximation reads: 
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For the one-loop corrections from the other sectors of the MSSM the logarithmic ap- 

In order 



proximation given in Ref. [12| has been used for FeynHiggsFast, see also Ref. [10 



to obtain the radiatively corrected Higgs boson masses and the mixing angle, the renor- 
malized Higgs boson self-energies have to be inserted into eq. (|6]), and the corresponding 
diagonalization has to be performed. 



We have also implemented the result of the MSSM contributions to Ap HS,|13"1. Here the 
corrections arising from t/b- loops up to 0(aa s ) have been taken into account. The result is 
valid for arbitrary parameters in the t- and 6-sector, also taking into account the mixing in 



the 6-sector which can have a non-negligible effect in the large tan {3 scenario |L3 



The two-loop result can be separated into the pure gluon-exchange contribution, which 
can be expressed by a very compact formula, allowing a very fast evaluation, and the pure 
gluino-exchange contribution, which is given by a rather lengthy expression. The latter 
correction goes to zero with increasing gluino mass and can thus be discarded for a heavy 
gluino. Concerning the implementation of the p -parameter into FeynHiggsFast, we have 
neglected the gluino exchange contribution. The p -parameter can be used as an additional 
constraint (besides the experimental bounds) on the squark masses. A value of Ap outside the 
experimentally preferred region of Ap SUSY m 1CT 3 \TA\ indicates experimentally disfavored 
t- and 6-masses. 



the complete leading two-loop 
see eq. (|) and 



For illustration of the quality of the compact approximation, we compare rrih and a e g 
with the results from the complete Feynman-diagrammatic calculation. The FD calculation 
contains the full diagrammatic one-loop contribution ||T5| 
corrections in 0(aa s ) HU], an d the two contributions beyond 0(aa, 
eq. (|I5D) without approximation. 

For the Standard Model parameters we use M z = 91.187 GeV, M w = 80.39 GeV, G F = 
1.16639 x 10 -5 GeV -2 , a s {m t ) = 0.1095, and m t = 175 GeV. In the numerical evaluation 
we have furthermore chosen the trilinear couplings in the scalar top and bottom sector to be 
Ab = A t . This fixes together with the choice of p (the Higgs mixing parameter) the mixing 
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Figure 1: nih as a function of rriq, calculated from the full formula and from the approxi- 
mation formula for M A = 500 GeV, m g = 500 GeV and tan (3 = 1.6 or 40. 




100 200 300 400 500 100 200 300 400 500 



M a [GeV] M A [GeV] 



Figure 2: m h as a function of M4, calculated from the full formula and from the approxi- 
mation formula for = 1000 GeV, m g = 500 GeV and tan/3 = 1.6 or 40. 
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Figure 3: The relative difference (sino; c ff — sino; e fj(approx))/ sma c s is shown as a function 
of Ma for three values of tan/5 in the no mixing and the maximal mixing scenario. The 
other parameters are \i = —100 GeV, M = rriq, m g = 500 GeV and A b = A t . 
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in the scalar bottom sector. The parameter M appearing in the plots is the SU(2) gaugino 
mass parameter, it enters as an independent parameter in the full result only, see Ref. 

In Refs. @,|10[ it has been shown that the lightest Higgs boson mass rrih as a function of 
Xt reaches a maximum at around \X t /rriq\ ~ 2. This case we refer to as 'maximal mixing'. 
A minimum of nth is reached for X t ~ 0, which we refer to as 'no mixing'. 

Fig. |l| displays the dependence of rrih on rrig for the cases of no mixing and maximal 
mixing, and we have set Ma = 500 GeV. For tan/5 we have chosen two typical values: 
tan/3 = 1.6 as a low and tan/3 = 40 as a typical high value. Very good agreement is found 
in the no-mixing scenario as well as in the maximal-mixing scenario, the deviations lie below 
2 GeV. 

The dependence on Ma is shown in Fig. ^|. The quality of the approximation is typically 
better than 1 GeV for the no-mixing case and better than 2 GeV for the maximal-mixing 
case. Only for very small (and experimentally already excluded) values of Ma a deviation 
of 5 GeV occurs. The peaks in the plot for tan/3 = 1.6 in the full result are due to the 
threshold Ma = 2m t in the one-loop contribution, originating from the top-loop diagram 
in the A self-energy. This peak does not occur in the approximation formula (where the 
momentum dependence of the A self-energy has been neglected) and can thus lead to a 
larger deviation around the threshold. 

In Fig. |3] we display the quality of the short analytical approximation formula for the 
effective mixing angle. The relative difference between a e g obtained from the full calcula- 
tion (i.e. calculating the renormalized Higgs self-energies S s (0) without any approximation) 
and the angle obtained with the help of the approximated Higgs self-energies, denoted as 
a e ff (approx) is shown as a function of Ma- We use three values of tan/3, tan/3 = 3,20,40. 
Sizable deviations occur only in the region 100 GeV < Ma < 150 GeV, especially for large 
tan /3. In this region of parameter space the values of rrih and ran are very close to each other. 
This results in a high sensitivity to small deviations in the Higgs boson self-energies entering 
the Higgs-boson mass matrix (^), (0). Otherwise the relative difference stays below 3%. 



3 The Fortran program FeynHiggsFast 

The complete program FeynHiggsFast consists of about 1300 lines Fortran code. The exe- 
cutable file fills about 65 KB disk space. The calculation for one set of parameters, including 
the Ap constraint, takes about 2 x 10~ 5 seconds on a Sigma station (Alpha processor, 600 
MHz processing speed, 512 MB RAM). The program can be obtained from the FeynHiggs 
home page: 

tittp : / /www-itp . physik . uni-karlsruhe . de/f eynhiggs . 

Here the code is available, together with a short instruction, information about bug fixes 
etc. 

FeynHiggsFast consists of a front-end (program FeynHiggsFast) and the main part 
where the evaluation is performed (starting with subroutine f eynhiggsf astsub). The 
front-end can be manipulated by the user at will, whereas the main part should not be 
changed. In this way FeynHiggsFast can be accommodated as a subroutine to existing 
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programs^ , thus providing an extremely fast and for many purposes sufficiently accurate 
evaluation for the masses and mixing angles in the MSSM Higgs sector. 

The input of FeynHiggsFast are the low energy SUSY parameters, listed in Tab. |I|. Con- 
cerning the t-sector, the user has the option to enter either the physical parameters, i.e. the 
masses and the mixing angle (to£ x , m~ t2 and sin 6^) or the unphysical parameters M^ L , Mi 
and X t . From these input parameters FeynHiggsFast calculates the masses and the mixing 
angle of the neutral CP-even Higgs sector, as well as the mass of the charged Higgs boson 
and the p -parameter. 



input for FeynHiggsFast 


MSSM expr. 


internal expr. in FeynHiggsFast 


tan (beta) 


tan/3 


ttb 


Msusy_top_L 


M h 


msusytl 


Msusy_top_R 




msusytr 


MtLR 


X t EE M t LR 


mtlr 


MSt2 




mst2 


delmst 


Am t - = m- t2 - m h 


delmst 


sin(theta_stop) 




stt 


MT 


m t 


ramt 


Mue 


p 


mmue 


MA 


M A 


mma 



Table 1: The meaning of the different MSSM variables which have to be entered into Feyn- 
HiggsFast. 



4 Conclusions 

FeynHiggsFast is a Fortran code for the calculation of masses and mixing angles in the Higgs 
sector of the MSSM. In addition it evaluates the SUSY corrections to the p -parameter, 
arising from the scalar top and bottom sector. Concerning the evaluation in the neutral 
CP-even Higgs sector, FeynHiggsFast is based on a simple analytic formula, derived from a 
more complete result obtained in the Feynman-diagrammatic approach. Beyond the one-loop 
level, FeynHiggsFast contains the dominant corrections in 0(aa s ) and further subdominant 
contributions. The accuracy of the approximation compared to the full result is better than 
2 GeV for most parts of the MSSM parameter space. 

2 This has been carried out, for example, for the program HDECAY |l6|,[l7]], into which FeynHiggsFast 
has been incorporated as a subroutine. 
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The program is available via the WWW page 
trttp : / /www-itp . physik . uni-karlsruhe . de/f eynhiggs . 

The code consists of a front-end and a subroutine. The front-end can be manipulated at 
the user's will. By accommodating the front-end, FeynHiggsFast can serve as a subroutine 
to existing programs, thus providing an extremely fast and reasonably accurate evaluation of 
the masses and mixing angles in the MSSM Higgs sector. These can then be used as inputs 
for further computations. 
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